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1. Introduction 

For the past decade there has been intense activity in the field of VT quantum mechanics 
[H 12]. A PT-symmetric Hamiltonian is said to have an unbroken VT symmetry if all 
of its eigenf unctions are also VT symmetric. A Hamiltonian having an unbroken VT 
symmetry is physically relevant because all of its eigenvalues are real and it generates 
unitary time evolution. Thus, such a Hamiltonian defines a conventional quantum- 
mechanical theory even though it may not be Dirac Hermitian. (A linear operator 
is Dirac Hermitian if it remains invariant under the combined operations of matrix 
transposition and complex conjugation.) One can regard such non- Hermitian quantum- 
mechanical systems as being complex extensions of conventional quantum systems. 

The interesting features of VT quantum mechanics have motivated many recent 
studies of VT classical mechanics. In particular, solutions to Hamilton's equations have 
been examined for various systems whose Hamiltonians are VT symmetric. For such 
systems the classical trajectories are typically complex [Sj HJ |5l El [71 El [U [TOl [TTl [121 II31 
[HI 1151 US]- These trajectories can lie in many-sheeted Riemann surfaces and often have 
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elaborate topological structure. When the VT symmetry of the quantum Hamiltonian 
is not broken, the real-energy trajectories of the corresponding classical Hamiltonian are 
found to be closed and periodic [HI fT3] . 

The purpose of this paper is to explore a new aspect of complex classical mechanics, 
namely, the complex extension of chaotic behavior. Specifically, we study two classical 
systems: the kicked rotor and the double pendulum. The kicked rotor is a paradigm for 
studying the dynamics of chaotic systems described by time-dependent Hamiltonians 
[T71 flSl [19] . The planar double pendulum is also a dynamical model whose classical 
motion is known to be chaotic [20]. The Hamiltonians for both of these dynamical 
systems are VT symmetric so long as the parameters K in the Hamiltonian for the 
kicked rotor ^ and g in the Hamiltonian for the double pendulum (12) are real. We 
use a variety of computational tools in order to derive the numerical results presented. 
The C programming language was used to implement of a fully symplectic three-stage 
Gauss-Legendre Runge-Kutta method for the simulation of the double pendulum, and 
standard functionality in Mathematica 6 was used in the study of the kicked-rotor. 

This paper is organized as follows: In Sec. [2] we define the kicked rotor and mention 
briefly the transition associated with the disappearance of KAM trajectories. In Sec. |3] 
we describe the planar double pendulum and describe the analogous transition that 
occurs for this dynamical system. We also reproduce the numerical work of Heyl 
concerning flip times. This work reveals fractal-like structure in the plane of initial 
conditions [21] . Then, in Sees. |4] and [s] we study the short- and long-time behaviors 
of the kicked rotor and the double pendulum in the complex domain, where in part 
our objective is to identify indicators for the transition to chaos. We also demonstrate 
that these two and very different dynamical systems exhibit remarkably similar features. 
Section |6] contains some concluding remarks. 



2. Kicked Rotor 



The Hamiltonian for the kicked rotor is [TTJ [181 ttS] 

H =^ + K COS e ^ 5{t-nT), (1) 

n=— oo 

where / is the moment of inertia of the rotor, p is its angular momentum, and 9 is the 
angular coordinate. As the rotor turns, it is subjected to a periodic impulse, which is 
applied at times t = 0, ±T, ±2T, .... The magnitude of the impulse is proportional 
to K, a constant having dimensions of angular momentum. This Hamiltonian is VT 
symmetric because it is symmetric separately under the operation of angular reflection 
V, where V : 6 ^ 27r ~6 and V : p —p, and the operation of time reversal T, where 
T : t ^ —t, T : p —p, and T leaves 6 invariant. (Note that angular reflection V is 
not the same as spacial reflection, which maps 9 ^ n — 9.) 
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Hamilton's equations of motion derived from ([T]) are 

? = ? and ^=Ksmey 5{t-nT). (2) 
dt I dt ^ ^ ^ ^ 

n=0 

These equations imply that the angular momentum p changes discontinuously at each 
kick, but remains constant between kicks. As a result, the angle 6 changes linearly with 
time if: between kicks and is continuous at each kick. 
It is customary to denote 

p„=p(nT + 0+) and 9^ = 9{nT + (3) 

Thus, Pn is the angular momentum and 6n is the angle variable immediately after the 
nth. kick. These variables satisfy the discretized version of ([2]): 

T 

9n+i =9n + jPn and Pn+1 = Pn + K smOn+i- (4) 

It is conventional to replace p„ and K by the dimensionless quantities 
T T 

jPn Pn and jK ^ K, (5) 
in terms of which we rewrite ^ in dimensionless form as 

On+i =On+Pn and pn+i = Pn + K siu 9n+l. (6) 

This system of difference equations, which depends on a single dimensionless real 
parameter K, is known as the standard map. It is straightforward to show that the 
standard map is area-preserving in p — 6 phase space. Note that the angular variable 
9n may be taken modulo 2%. It then follows from the first equation in ^ that Pn may 
also be taken modulo 2tc. Thus, (|6| maps the two dimensional torus onto itself. 

The behavior of the standard map ^ is elaborate and has been studied extensively 
[nillHlIiniESlEalEllESlEH]. For small K the motion in phase space is bounded and 
chaotic in some regions. As K increases, KAM trajectories disappear. At the critical 
value Kc = 0.9716 . . . only the KAM trajectories with golden-mean winding number and 
with inverse-golden-mean winding number remain, and the motion in phase space is still 
confined. For K > Kc the last bounding trajectory is destroyed and global diffusion in 
phase space ensues. The critical behavior near Kc has been studied intensively [231126] . 

Figure [T] illustrates the transition from subcritical to supercritical K for the kicked 
rotor. In this figure we display four sets of superpositions of phase planes, each consisting 
of eleven randomly chosen initial conditions 6o,po- For each set of initial conditions we 
allow the time variable n to range from 1 to several thousand. The values of K for these 
four plots are 0.40, 0.97, 2.0, and 4.0. 

In this paper we continue the classical dynamics described by the standard map ^ 
into complex phase spacejf] Our objective here is to generalize ^ into complex phase 
space and thereby gain a better understanding of the critical behavior near Kc. To 

^ The idea to study chaotic systems in complex phase space was introduced several years ago in 
Ref. [27] . The motivation in these papers was to study the effects of classical chaos on semiclassical 
tunneling. In the instanton calculus one must deal with a complex configuration space. 




Figure 1. Four phase-plane views of the kicked rotor. In each of the figures we display 
the superposition of the discrete phase-plane trajectories 9m Pn for eleven randomly 
chosen initial conditions. The time variable n ranges from 1 to several thousand. The 
values of K for the four plots are 0.40, 0.97, 2.0, and 4.0. The KAM surfaces separate 
Poincare islands. Observe that as K increases, the KAM surfaces gradually disappear 
and the trajectories diffuse into the phase plane. 



accomplish this we are motivated to extend the analysis of Refs. [31 SI El El [71 El [9l UHl 
HH [El [131 mi [m US] to time-dependent systems. Thus, we treat pn, On, and sometimes 
K as complex variables, which we separate into real and imaginary parts as 

Pn = rn + iSn, 6n = On + iPn, K = L + iM. (7) 

Substituting ([T]) in (|6]), we obtain the complexified standard map 

r„+i =rn + L sin cosh (3n+i - M cos sinh 

COS sinh M sin cosh (8) 

In Sees. |4]and|5]we display and discuss the results of our numerical studies of (|8]). 
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Figure 2. Configuration of the double pendulum. The double pendulum consists of 
two massless rods, each having a massive bob at the end. The second rod hangs from 
the end of the first rod. The two rods are constrained to swing in a plane and are 
acted on by a homogeneous gravitational field of strength g. 



3. Double Pendulum 

As shown in Fig. |2| a planar double pendulum consists of a massless rod of length ii 
with a bob of mass mi at the lower end from which hangs a second massless rod of 
length £2 with a second bob of mass 1712 at the lower end. This compound pendulum 
swings in a homogeneous gravitational field g, and its motion is constrained to a plane. 

In this paper we take both bobs to have unit mass and both rods to have unit 
length. The coordinates of the bobs in terms of the angles from the vertical are 

xi = sin^i, 

yi = - cos 6*1, 

X2 = sin 9i + sin 62, 

1/2 = — cos — cos 6*2- (9) 

Therefore, the potential and kinetic energies of the double pendulum are 

= -^cos02 -2^cos^i, and T = Ol + + cos{9i ~ 62) . (10) 

From these one can form the Lagrangian L = T — V, and then construct the Hamiltonian 
for the system by a Legendre transform. We obtain 

rr pI + - 2piP2 COs(gi - 62) n n 

H = p — TT, ^ 1 q COSO2 — 2q cosOi. 11 

2 [sin^Oi - 62) + I] y 2 y 1 k ) 

This Hamiltonian is PT-symmetric because it is symmetric separately under the 
operation of angular reflection V, where V : 9i 2 ^ 2-71 — 61^2 and V : pi,2 — >■ — Pi,2, and 
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the operation of time reversal T, where T : t —>■ —t, T : pi^2 —Pi,2, and T leaves 61^2 
invariant. 

Hamilton's equations are then 

oOi sm (^1 -62) + I 

[pI + 2pl - 2pip2 cos(^i - ^^2)] sin [2{9i - ^^2)] 



+ 



2[sin2(^i-^2) + l]' 



dH ■ n , PiP2sm 
P2= - = -g sm 02 + Z772 



01 



802 sin' (^1 - ^2) + 1 

[pj + 2pl - 2pip2 cos {01 - 02)] sin [2 (^i - ^2)] 

2[sin2(^i-^2) + l]' 
dH pi - p2 cos {01 - 02) 



02=-£r= ^i~2.r • (12) 



dpx sin^ - ^2) + 1 
dH_ ^ 2p2 - Pi cos (gi - 02 ) 
dp2 ~ sin' - + 1 
Note that this system conserves energy, unlike the kicked rotor whose Hamiltonian ([T]) 
is time-dependent. 

A beautiful and convincing numerical demonstration that the motion of the double 
pendulum is complicated and elaborate was given by Heyl [21]. In his work Heyl 
calculates for a given initial condition the time required for either pendulum to exhibit a 
flip\ that is, for either 0i or 02 to exceed the value vr. This calculation is then performed 
for the limited set of initial conditions for which pi(0) = and ^2(0) = 0, and the initial 
values of 0i and 02 both range from to 27r. Each pixel in the initial 6'i, 02 plane is then 
colored according to the length of the flip time. 

We have applied Heyl's approach to the double pendulum in ( 12 ) and have used a 
Gauss-Legendre Runge-Kutta method, which is known to be fully symplectic [281 129] . 
The results of this calculation are given in Fig. |3} This figure is composed of a 600 x 600 
grid of pixels, where each pixel represents the initial condition [6'i(0), 6'2(0)]. If neither 
pendulum flips within 100 time units, then the pixel is assigned the color black (the 
convex-lens-shaped region in the center of the figure). If either pendulum flips in a 
short time, the pixel is colored dark gray, with longer flip times being indicated by 
lighter shades of gray. Notice the fractal-like structure throughout the diagram. The 
appearance of this complicated structure demonstrates that even though the double 
pendulum has only two degrees of freedom, it exhibits rich and nontrivial dynamics. 

In analogy with the kicked rotor, there is a transition in the behavior of the double 
pendulum in which KAM surfaces disappear as a dimensionless parameter increases 
beyond a critical value. This parameter, which measures the strength of the gravitational 
field relative to the total energy, is defined as [20] 

7 = m,gii/E. (13) 

The transition occurs near 7 = 0.1 and Ref. [20] shows that at the transition the 
last surviving KAM surface is the one with winding number being equal to the golden 
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Figure 3. Flip times for initial conditions — tt < Oi{0) < tt, — tt < 02(0) < tt, pi(0) = 0, 
and P2(0) = 0. By a flip we mean that the angular position of either bob exceeds tt. 
This figure contains a 600 x 600 grid of pixels. The color of each pixel characterizes the 
behavior of the double pendulum that arises from the initial condition [0i(O), ^2(0)]. If 
neither pendulum flips within 100 time units, then the pixel is black (the convex-lens- 
shaped region in the center of the figure). If one of the pendula flips in a short time, 
the pixel is colored dark gray. Longer flip times are indicated by lighter shades of gray. 
The fractal-like structure throughout the diagram reveals the nontrivial dynamics of 
the double pendulum. Because of parity symmetry this figure is symmetric under 
the combined reflections 6*12 —6*1, 2- Note that different pixels are associated with 
different energies and that there is an elliptical region at the center of the figure in 
which flips are forbidden by energy considerations. 

mean. In Fig. |4]we plot the Poincare sections generated from 25 randomly chosen initial 
conditions for four different values of 7. The plot displays points in the 61, pi plane when 
62 = and simultaneously ^2 > 0. A KAM surface is visible when 7 = 0.05, which is 
below the critical value. At 7 = 0.1, which is near the critical value, the KAM surface 
disappears. For the other two values of 7, which are significantly greater than the critical 
value, the distribution of points in the plot becomes diffuse in a manner analogous to 
the behavior displayed in Fig. [T]for the kicked rotor in the K > Kc regime. 

4. Short-time behavior 

Having reviewed some of the well-known properties of the kicked rotor and the double 
pendulum, we now proceed to examine the behavior of the solutions to the kicked-rotor 
and double-pendulum equations of motion in the complex domain. To do so, we do 
not change the form of the equations of motion, but rather we take complex initial 
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Figure 4. Poincare plots for the double pendulum for four values of 7. Below the 
critical value, which is near 7 = 0.1, the final KAM surface can still be seen in the 
second plot, but this surface evaporates as 7 increases past the critical value and the 
points in the plot spread. This plot is analogous to that in Fig. [T] for the kicked rotor. 



conditions and in some cases we allow the parameter in ([6]) for the kicked rotor and 
the parameter in (12) for the double pendulum to take on complex values. 

In this section we investigate the behavior of these dynamical systems for short 
times; that is, for up to 1000 time steps. For the kicked rotor, let us see what happens 
if we take the initial momentum to be real, po = 0, but take the initial angle to have 
a small imaginary component, 6q = 1 + O.OOOli. In Fig. |5] we plot the points On in 
the complex-^ plane for n = 0, 1, 2, . . . , 1000 for a range of real values of K around 
the critical point Kc = 0.9716 .... While it is difficult to see the subtle change from 
subcritical to supercritical behavior in real plots like those in Fig. [T| a qualitative change 
in the complex behavior is quite evident in Fig. [s] pO]. Below the points tend to 
occupy a two-dimensional region in the complex plane with some fine structure in it, but 
as K increases above points tend to coalesce along well separated one-dimensional 
curves. 

In analogy with Fig. |5} we plot in Fig. [6] a trajectory for the double pendulum in 
the Re^i, Imp2 plane as a function of t for < t < 224. We take two values of 7, one 
subcritical and one supercritical and use the shghtly complex initial conditions 9i = 3.1, 
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Figure 5. Behavior of the solution to the kicked rotor in the coniplex-6' plane for a 
range of real values of K near Kc ~ 0.9716 .... We take as an initial condition po = 
and 6q — 1 + O.OOOli and allow the system to evolve for n — 1000 time steps. For each 
value of n we plot the complex value of 0„ as a point in the complex- plane. The 
plot corresponding to the critical value is highlighted. Note that there is a qualitative 
change in the nature of these complex plots as K passes through its critical value. 
Specifically, the points making up the plot become less uniform and more stratified. 
These changes in behavior are easier to observe than those in Fig. [T] 



62 = 3.1 + O.OOOli, pi = 1.283, and p2 = 1.283 in which the numbers are chosen at 
random. Observe that when 7 is above the critical value the trajectory appears to be 
confined to distinct narrow bands, somewhat similar to the stratified structures that 
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occur for K > Kc in Fig. |5] for the kicked rotor. However, when 7 is below the critical 
value, the trajectory spreads and is similar to the more diffuse behavior in Fig. [5] when 
K < Kc- It would be interesting to understand this stratification analyticallty and also 
to understand its relation to KAM theory. 

An even more dramatic way to observe the transition from the subcritical to the 
supercritical regions of the kicked rotor is to construct plots like that in Fig. [3] for the 
double pendulum. We take as an initial condition po = and take a range of complex 
initial values for Oq: — tt < Re^o ^ and — vr < Im^^o ^ tt. For each initial condition we 
then allow the kicked rotor to evolve up to a maximum of 400 time steps and determine 
the time step at which the real part of p„ becomes infinite, if it does actually become 
infinite. (Here, by infinite we mean that the numerical value of Kepn exceeds 10^°^, 
which is the largest number that may be represented in double precision arithmetic.) 
We then perform this calculation for each pixel on a 628 x 628 grid representing the 
complex-^o plane. We assign a color to each pixel corresponding to the time at which 
Pn becomes infinite: White indicates that pn does not become infinite within 400 time 
steps, and darker shades indicate that p„ becomes infinite after shorter times. 

In Fig. [7] we display the results of this calculation for K = 0.01, 0.1, 0.2, and 0.6, 
and in Fig. [8] we display the results for this calculation for K = 0.9, 1.1, 2.0, and 5.0. 
Note that all these figures exhibit a complicated dendritic and fractal-like structure. 
A number of qualitative changes occur as K increases past its critical value. One 
obvious change is that the dendritic landscape becomes smoother and more rounded as 
K increases. A less obvious change is that the regions in which Rep„ does not diverge 
for n < 400 become connected when K exceeds Kc. 

Instead of constraining K to be real, we can, of course, take K to be complex. In 
Fig. |9]we take K = 0.6i. Note that the fractal structure in Figs. [7] and |8] is preserved, 
but it is distorted and loses its left-right symmetry. 

Rather than requiring that K pass its critical value on the real-i^ axis, it is possible 
to go from a subcritical real value to a supercritical real value via a path in the complex- 



K plane, as in Fig. 10 The pictures making up this figure are constructed from 
values of K that lie on a semicircle of radius 0.6 and are centered at the critical value 
Kc = 0.9716 .... In these figures we observe fractal structures like those in Figs. [^{9} 
but they are slightly distorted. However, there is a significant difference in that there 
is mottling (replacement of large patches of solid shading by a speckled pattern) in 
the graphs where K is complex; this mottling is absent when K is pure real or pure 
imaginary. (Note that when K is complex, VT symmetry is broken if the T operator 
is antilinear, that is, it changes the sign of zj§]) 

§ When K and g become pure imaginary the system becomes invariant under combined VT reflection. 
However, now is a spatial reflection, V : 9 ^ 6 + tt,so that both cos 6 and sin 6*, and thus the cartesian 
coordinates, change sign. The sign of the angular momentum now remains unchanged under parity 
reflection. This explains the symmetry of the plots when K and g are pure imaginary (see Fig. |9] and the 



lower-right plot in Fig. 11 respectively). This change of symmetry of the system as its couplings vary 
in complex parameter space is not unusual. For example, at a generic point in coupling space for the 
three-dimensional anisotropic harmonic oscillator, the only symmetry is parity. However, when any two 
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Figure 6. Trajectory for the double pendulum that results from the slightly complex 
initial conditions di — 3, 02 = i + O.OOOli, pi ~ p2 — 2. The parameter 7 is varied 
by changing the values of g and li in relation (13 1 and fixing mi — E — I. This plot 



displays the trajectory in the Re^i, Imp2 plane as a function of t for < i < 1000. 
The trajectory is confined to narrow bands when 7 is above the critical value and is 
more diffuse when it is below the critical value. In this sense, this figure has features 
that are similar to those in Fig. [5] for the kicked rotor. 



couplings coincide and are different from the third, the reflection symmetry is enhanced and becomes 
a continuous symmetry, namely, an 0(2) symmetry around the third axis. (There remains parity- 



Chaotic systems in complex phase space 
K=0.01 



K=0.1 



12 




Figure 7. Behavior of the solution to the kicked rotor in the complex-6' plane for four 
real values of K and a range of complex initial conditions: po = 0, < Re^o ^ 27r, 
and < Im^o < 27r. We allow the system to evolve for at most n — 400 time steps. 
We assign a color to each value of 9o according to the time at which Rep„ becomes 
infinite (if it does so). The graph of initial values of 9 clearly has fractal structure. 
Distinctive changes occur in the nature of these complex plots as K increases. 

We have performed a closely related study for the double pendulum: We allow 
g to be complex and repeat the numerical analysis of the short-time behavior that 
we used to produce Fig. [3j We find that as the imaginary part of g increases, 
the fractal-like structure that we see in Fig. |3] gradually moves outward towards the 

time reflection symmetry in the third direction.) When all three couplings coincide, the symmetry is 
enhanced further and becomes a full 0(3) symmetry. 
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Figure 8. Same as in Fig. [7]with four higher values of K. There are many quahtative 
changes that occur as K increases. For example, the dendritic landscape that occurs 
for smaller values of K becomes smoother and more rounded as K increases. A subtle 
but important change is that the white regions, the regions in which Rep„ does not 
diverge for n < 400, become connected when K exceeds Kc- 



corners of the figure. Correspondingly, the boundaries between different colored regions 
become smoother. To demonstrate this effect, we choose four different complex values, 
^ = 0.1 + O.OOSi, 0.1 + OMi, 0.1 + O.li, O.li, and plot the results in Fig. In 
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Figure 9. Same as in Figs. [7] and [8] except that K = 0.6i. Observe that while many 
features of the fractal structure in Figs. |7] and [8] are preserved, they are distorted and 
the left-right symmetry is destroyed. 



5. Long-time behavior 

In this section we study the long-time behavior of the kicked rotor and the double 
pendulum in complex phase space and we find that they share many qualitative features 
in this regime as well. By long-time we mean roughly 10^ to 10^ time steps or time units 
rather than the 10^ time steps taken in Sec. |4} We find that the solutions to the 
dynamical equations for these systems exhibit characteristic behaviors at different time 
scales. On a long-time scale, which is determined by the imaginary part of the initial 
value of an angle, the solutions tend to ring; that is, the envelope of the solution grows 
and decays to zero with gradually changing periods. On a short-time scale the solution 



exhibits a distinct and clearly identifiable rapid oscillation, as we can see in Fig. 12 
We have chosen here to use the language of multiple-scale perturbation theory 
to describe this oscillatory behavior. However, for both the kicked rotor and the double 
pendulum the unperturbed equations are not linear, and thus the usual techniques of 
multiple-scale perturbation theory cannot be applied directly in these cases. 

Let us first examine the kicked rotor. As an initial condition we choose po = and 
6*0 = 1 + ie. In Fig. 12 we take K = 0.6 and e = 10~^ and we plot Re6'„, Im6'„, Repn, 



and Impn for < n < 260 000. Note that while Re6'„ and Rep„ oscillate within almost 
constant boundaries, Im6'„ and Imp„ appear to ring with a period of order 1/e. To 
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Figure 10. Same as in Fig. [7] except that the pictures making up this figure are 
constructed from values of K lying on a semicircle of radius 0.6 and centered at the 
critical value Kc — 0.9716 ... .In these figures we observe fractal structures similar to 
those in Figs. |7]|9] but slightly distorted. An important difference between this figure 
and Figs. 7][9 is that there is mottling (speckling) in those graphs where K is complex. 



verify this dependence on e we take e = 10~^, which is ten times larger, and we do not 
change the other initial conditions or the value of K. The result for Im^^ is given in 
Fig. [131 where we see that the period of the ringing is roughly ten times shorter than the 



period in Fig. 12 In both of these figures the ringing eventually comes to an abrupt end, 
at which point the iteration diverges and the amplitude of oscillation becomes infinite; 
this happens after about 2| rings in Fig. 



12 



and after about eleven rings in Fig. 



13 



While the inverse of the imaginary part of appears to set the scale of the ringing 
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Figure 11. Short-time behavior of the double pendukim with complex g. Same 
as Fig. [3] but with g = 0.1 + O.OOSi, 0.1 + O.OH, 0.1 + O.lz, O.li. Note that as Im^ 
increases, the fractal-like structure that we see in Fig. [3] moves outward and towards 
the corners, and the boundaries between differently colored regions become smooth. 
In contrast with Fig. |3] there are no energetically-forbidden-flip regions because there 
exist complex pathways from any pixel to a flipped configuration. 



period, we have found that the length of the ringing period is also sensitive to the value 
of K. In Fig. 14 we plot Impn for four values of K: 0.5, 0.53, 0.55, and 1.0. For each of 
these values of K we plot Imp„ until it diverges. The initial conditions for each graph 



in Fig. 14 are the same as those in Fig. 12 



The double pendulum exhibits a long-time ringing behavior that almost exactly 
parallels that of the kicked rotor. We plot the long-time behavior of Im^i for = 1 
and initial conditions pi{0) = ^2(0) = 0, ^^i(0) = 1, and 6*2(0) = 10~^i in Fig. 15 and for 
^2(0) = 2 X 10~^2 in Fig. 16 Note that, like the kicked rotor, the long-time-scale ringing 
periods are determined by the imaginary part of the initial value of an angle; here, the 
period is proportional to 1/ [Im^2( 
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6. Concluding remarks 

Apart from making the obvious remark that the two nonlinear systems studied in this 
paper exhibit very similar short-time and long-time dynamical behaviors, this work 
indicates that studying the dynamics of classical chaotic systems in complex phase space 
may help us to understand the onset of chaos. For example, for the case of the kicked 
rotor, we observe in Fig. |5]a change in the complex behavior as K increases past K^. Of 
course, the results reported here are empirical, but they clearly underscore the need for a 
deeper analytical understanding of these models. For example, an important unanswered 
question is, what is the analog of the KAM theorem in complex phase space? 

Finally, we remark that the kicked rotor is one of the rare time-dependent systems 
whose quantum dynamics may be studied in some detail. Indeed, the kicked rotor is a 
paradigm for studying quantum chaos. It might be particularly useful to explore the 
PT-deformed analog of the work of Fishman et al. [19] because (i) this would be a 
nontrivial extension of VT quantum mechanics to time-dependent systems, and (ii) it 
may provide a way to define and understand PT-symmetric quantum chaos. 
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Figure 12. Long-time behavior of the kicked rotor with initial conditions pq — 
and 9o = 1 + lO^^i and K = 0.6. Note that while the real parts of p„ and On 
oscillate between almost constant boundaries, the imaginary parts of Pn and 9n exhibit 
a synchronized ringing behavior whose period is of order 1/ (Im^o)- In addition to 
the long-time ringing there is a short-time oscillation that becomes most pronounced 
when the amplitude of the ringing is at a maximum. After about 2^ rings the solution 
abruptly diverges and ceases to exist. 
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Figure 13. Same as in Fig. 12 except that the imaginary part of Oq is ten times 
smaller: = 1 + lO^'^z and only the imaginary part of is displayed. In this figure 
the period of the ringing is ten times shorter and there are about eleven rings before 
the solution destabilizes and ceases to exist. 
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Figure 14. Long-time behavior of Inip„ for the kicked rotor for four different values 
of K. The initial conditions are the same as in Fig. [12] Observe that the period of 
ringing is quite sensitive to the value of K. 
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Figure 15. Long-time behavior of the double pendulum with g = I. The initial 
conditions used for this plot are pi{0) = P2(0) = 0, Oi{0) = 1, and ^2(0) = lO^^i. 
Like the kicked rotor, the imaginary part of an angle exhibits a ringing behavior whose 
characteristic period is of order 1/ (Im02(O)). The plot terminates when the solution 
to the equations of motion abruptly diverges. Similar to the behavior of the kicked 
rotor, there is also a short-time oscillation, but unlike the kicked rotor, the positive 
and negative peaks are out of phase with one another. 
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Figure 16. Same as Fig. [Ts] but with ^2(0) = 2 x IQ-H. Note that doubling the 
imaginary part of 6*2(0) has the effect of roughly halving the ringing period. 



